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Abstract. A new and computationally viable full quantum version of line shape 
theory is obtained in terms of a mixed Weyl symbol calculus. The basic ingredient in 
the collision-broadened line shape theory is the time dependent dipole autocorrelation 
function of the radiator-perturber system. The observed spectral intensity is the 
Fourier transform of this correlation function. A modified form of the Wigner- 
Weyl isomorphism between quantum operators and phase space functions (Weyl 
symbols) is introduced in order to describe the quantum structure of this system. 
This modification uses a partial Wigner transform in which the radiator-perturber 
relative motion degrees of freedom are transformed into a phase space dependence, 
while operators associated with the internal molecular degrees of freedom are kept 
in their original Hilbert space form. The result of this partial Wigner transform 
is called a mixed Weyl symbol. The star product, Moyal bracket and asymptotic 
expansions native to the mixed Weyl symbol calculus are determined. The correlation 
function is represented as the phase space integral of the product of two mixed symbols: 
one corresponding to the initial configuration of the system, the other being its time 
evolving dynamical value. There are, in this approach, two semiclassical expansions — 
one associated with the perturber scattering process, the other with the mixed symbol 
star product. These approximations are used in combination to obtain representations 
of the autocorrelation that are sufficiently simple to allow numerical calculation. The 
leading O(fi ) approximation recovers the standard classical path approximation for 
line shapes. The higher order Oih 1 ) corrections arise from the noncommutative nature 
of the star product. 
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1. Introduction 

Collision-broadened spectral line shapes carry important information on the time- 
dependent dynamical and collisional processes occurring in a radiating medium. For 
a system consisting of radiators (light emitting or absorbing species) and perturbers 
(atoms, molecules, ions or electrons), the basic ingredient in spectral line shape theory 
is the Fourier transform of the time autocorrelation function of the radiator dipole. The 
physical processes here involve the dynamical evolution of the radiator as it undergoes 
multiple collisions with a large number of perturbers. The interaction between the time 
dependent radiator dipole moment and the electromagnetic field induces absorption 
and emission. In this paper we develop a phase space based representation of quantum 
mechanics that is suitable for determining this dynamics. Within this formalism there 
are two different semiclassical approximations. These are employed in tandem to obtain 
representations of the autocorrelation function that are sufficiently simple to allow 
numerical calculation. 

The method developed here is suggested by the Moyal description [||] of quantum 
mechanics, which is based on the Wigner-Weyl isomorphism || between Hilbert space 
operators and functions (Weyl symbols) on classical phase space. In this isomorphism |J 
the canonical (qi,pi) operators are transformed into linear phase space functions (qi,Pi). 
The line shape problem presents one with two distinguished degrees of freedom: the 
internal molecular coordinates and the relative radiator-perturber separation variables. 
In many circumstances (large mass, large impact parameter, high relative velocity, weak 
pairwise interaction) one expects the relative motion to be nearly classical. On the 
other hand, evolutions of the subsystems associated with the internal coordinates are, 
generally, far from classical. This circumstance is exploited by the introduction of a 
partial Wigner transform that converts the relative motion degrees of freedom into 
phase space variables, but keeps operators associated with the internal coordinates in 
their original Hilbert space form. The result of this transformation we call a mixed Weyl 
symbol. 

Our approach to a quantum line shape theory is strongly influenced by the well 
established 'classical path approximation', which is used in many areas of atomic and 
molecular physics. It assumes that the perturbers move as classical particles, that is, 
along definite paths determined by the radiator-perturber interaction. The traditional 
justification for this approximation arises from the notion that one may consider the 
motion of perturbers in terms of packets of translational wave functions, following 
the laws of classical mechanics. Anderson^] in setting up his pioneering line shape 
theory argued, that since the typical distances of closest approach are about 5A, an 
uncertainty in position of lA leads to no great ambiguity in the magnitude or the type 
of intermolecular forces involved in the collision. He invoked the uncertainty principle to 
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conclude that the corresponding uncertainty in velocity is only a few percent for most 
molecular pairs. The approximation evidently breaks down || when the de Broglie 
wavelength of the perturber is comparable to or larger than a characteristic distance 
over which the intermolecular potential varies appreciably. An alternative way || of 
expressing this condition is that for the profile to be adequately treated with a classical 
path, the collisions of significance must have large angular momentum. Such a condition 
is reminiscent of that required to calculate scattering cross sections classically 0. Early 
derivations of the approximation employed heuristic arguments which relied on physical 
insight. The subsequent work of Baranger|| and Smith et al f| unified the derivation 
of the classical path approximation, and in particular employed statistical mechanics to 
justify the representation of the correlation function as an integral over the phase space 
variables of the perturber. The phase space average procedure || remains a cornerstone 
of the theory to the present. 

The system we shall consider is a dilute gas in which each of the radiator and 
perturber subsystems is in equilibrium and in which the binary collision approximation 
is valid. The binary collision approximation, which implies that strong collisions between 
radiator and perturber are well-separated in time, is a central assumption ]10 in most 
line broadening theories. The radiator density is sufficiently low that the radiator- 
radiator interaction may be ignored. The linewidth is dominated by collision broadening 
effects; Doppler and lifetime broadening are omitted from the analysis. The heavy 
radiator approximation is thereby imposed. The radiator and perturber subsystems are 
statistically independent; that is, the state of the perturber does not depend explicitly 
on the state of the radiator and vice versa. This assumption is frequently termed the 
'lack of back reaction' in the density matrix. Consequently the density matrix may be 
factored as p = p^p^ p \ where p( R > and p( p > depend only on the radiator and perturber 
variables, respectively 0] . Finally, since interest lies primarily in the effect of the field- 
radiator interaction on the radiator, the electromagnetic field is treated classically. For 
an overview of collision-broadening, the reader is referred to the review || of Allard 
and Kielkopf. 

A traditional starting point for the development of line shape theory is the expression 
for the power gained or lost from the radiation field to the molecular many body system 
in making electric dipole transitions from all initial states i to all possible final states 
/. The power produced by a single radiator interacting with the perturbers, is given 
by YtiffruifPiPif where is the Fermi Golden Rule probability per unit time for a 
transition between states i and / having an energy difference E{ — Ef = huoif. The 
weight factor p^ is the initial state density. The frequency content of this expression is 
known as the line shape or spectral profile, I(uj). If C^it) is the iV-perturber dipole- 
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dipole autocorrelation function, the spectral profile is the Fourier integral 

1 r°° 

I (to) = -Re / dte luJt C N (t). (1.1) 

7T JO 

Throughout the paper we employ the version of the binary collision approximation 
H 0, which expresses Cjsr(t) via the one perturber autocorrelation function C(t), cf. 



(OD, namely C N (t) = C{t) N . 

In this setting, the Hamiltonian required in the one perturber autocorrelation 
function is a sum of three contributions 

H = H 1 + H 2 + H 12 . (1.2a) 

The Hamiltonian Hi determines the (radiator) molecular structure, and the pair of 
Hamiltonians H 2 + Hi 2 generate different parts of the radiator-perturber dynamics. 
Let the d\ dimensional internal radiator coordinates be the Cartesian variables Q = 
(Qi, . . . , Qd-i) ^ R-q 1 an d the radiator-perturber coordinates given by q = (gi, . . . , q^ 2 ) G 
Rf. The state spaces over these two independent coordinate manifolds are Hi = 
L 2 (R dl ) and H 2 = L 2 {R d2 ). The full system Hilbert space is H = Hi <g> H 2 = 
L 2 (R dl xR d2 ). For reasons of clarity of presentation and notational convenience, we 
assume that the perturber may be treated as a point particle. This means that q is 
a three dimensional vector (d 2 = 3). The mixed Weyl symbol calculus and line shape 
theory developed here are valid for all values of d\ and d 2 . If desired, the one perturber 
Hamiltonian ( |1.2o| ) may be extended to a many perturber system. 

The Hamiltonians Hi and H 2 are simple in the sense that they are tensor products of 
operators on Hi and H 2 , namely Hi = hi ®I 2 and H 2 = Ii®h 2 , where Ii is the identity 
on H{. The Hamiltonian hi and its associated eigenvalue problem hi \Q n ) = E n \$ n ) 
determine the energy spectrum and wave functions of the radiator. The wave function 
|<3>i) is the molecular groundstate. If h 2 is the perturber kinetic energy and v 2 is a Q- 
independent (isotropic) part of the intermolecular interaction energy, then h 2 = h 2 fl + v 2 . 
The operator Hi 2 includes all the anisotropic parts of the radiator-perturber potential. 
It depends on the relative orientation of the molecular axis and the vector q. This 
operator is a function of both q and Q. For example in the case where the radiator can 
be treated as a rigid rotor (d x = 3), then if 12 has the Legendre polynomial expansion 

oo 

Hi 2 = Y,vS(\q\)m-Q)- (1-26) 

1=0 

The multipole potentials V^CM) are phenomenologically known functions. The 
presence of the term in ([1.2 fl) allows one to include portions of the isotropic potential 



in H 12 . 

The dipole autocorrelation function [|, [^, |1^] generated by the dynamics of H is 

C(t) = Tr w t)^ 3 U{H- t)p) . (1.3) 
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In the formula above, U(H;t) = exp(—iHt/h) denotes Schrodinger picture evolution 
and p,j is the jth Cartesian component of the radiator electric dipole. Tensors with 
repeated indices such as j G (1, 2, 3) are summed over. The initial (t = 0) thermal state 
of the system is specified by the density matrix p = e _/3 ^/Tre -/3H , where f3 = (&T) _1 . 
The statistical structure of the radiator-perturber system is normally assumed to be a 
stationary ergodic ensemble with temperature T. In this case the canonical ensemble 
average on the right-hand side of ( |1.3| ) is equivalent to the long time average over many 
radiator-perturber collisions. An attractive feature of the correlation function C(t) is 
that the effects of statistics (via p) and dynamics (via U(H;t)) are clearly separated. 

Our approach to computing C(t) is to generalize the existing Moyal quantum 
mechanics so that it can predict the H dynamics in the full Hilbert space Ti. In order 
to see the possible relevance of the Wigner transform methods to the classical path 
approximation, consider in isolation the one body (perturber) problem generated by /12 
on 7i 2 = L 2 (R d2 ). Adapted to this setting, the key elements || of the Moyal theory 
are as follows. Operators A on 7i 2 m ay be represented as functions (Weyl symbols) on 
phase space. The relevant phase space, with variables z = (q,p), is that induced by the 
manifold R^ 2 , namely T 2 * = T*(R^ 2 ) ~ R^ 2 x R^ 2 . We label the Wigner transform |]| 
map a, i.e. aA(z) = A w (z), where 

A w (z)= I dxe- ip - x/n (q + \x\A\q-\x). (1.4) 

The Fourier transform nature of the Wigner correspondence (|1.4j ) ensures that it has an 
inverse, cr _1 . 

The Heisenberg picture evolution generated by h 2 has two equivalent realizations. 
In Hilbert space 7^2, one has A(t) = T(h,2;t)A = U(h,2,tyAU(h2',t). In symbol space 
the equivalent of T(h 2 ;t) is denoted by r(h2',t) = aT{h2;t)a^ 1 . This latter operator 
maps symbols to symbols. In detail, if A w (z) is the Weyl symbol of A at t = and 
v4(t|z) = (A(t)) w (z) is the corresponding dynamical value at time t, then 

r(h 2 ;t)A w (z) = A(t\z). (1.5a) 

Acting on suitable symbols, r(h 2 ; t) defines a one parameter group with the same 
structure one finds with r(/i 2 ;t), namely, r(h 2 ;ti + t 2 ) = r(h 2 ;t2)r(h 2 ;ti). The 
equation of motion for -A(t|z) follows from taking the Wigner transform of the Heisenberg 
equation for A(t), 

jA(t\z) = {A(t),h2}u{z). (1.56) 

Here {■, -}m is the Moyal bracket, cf. ( ETTTJ ), of the symbol pair A{t) and h 2 = (^2)w 

Given the solution of ( |1.56p , expectation values in 7i 2 are realized as integrals over 
T 2 *. Let ip G 7^2 be a unit normalized initial state defining a density matrix p = \%l))(%p\. 
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In terms of the Wigner distribution w^(z) = (h)~ d2 (p) w (z), (h 



(A(t))^ = Tr H2 pA{t) 



i'1 



dzw^z) r(h 2 ; t)A^(z) 



'1.6a) 



In the circumstances where h 2 is ^-independent, 
asymptotic expansion 



r(h 2 ;t) admits a small h 



-2n 



r(h 2 ;t)A w = 77TTll {2n) (h 2 -t)A v 



n 2 



l { »\h 2 ;t) + -^\h 2] t) 



A w + 0(h 4 



[1.6b) 



n=0 M 

Like r(h 2 ;t), the quantities ^ 2n \h 2 ; t) are maps on symbols. In reference explicit 
formulas for ^ 2n \h 2 ; t) were derived as a result of a connected graph representation for 
r{h 2 ; t). The leading term "f^(h 2 ; t) (as has long been known [|14], [15], |l6l| ) is determined 
by the h 2 generated classical flow. The higher order operator coefficients, beginning 
with j^(h 2 ;t), have the form of partial differential operators [cf. (|3.7c )1 acting on 



symbols. Combining expansion ( |1.6£j ) and ( 1.6oj) yields a semiclassical expansion for 



the expectation value (A(t))$. With suitable restrictions on A w , rigorous error bound 
estimates are available for the asymptotic expansion (OHj). The early work of Antonets 
HU verified, for finite time displacements t G [0,T], that r(h 2 ,t)A w — > 7^ (h 2 ;t)A w in 
an appropriate norm as h —>■ 0. More recently new proofs []T7|] have been constructed 
that obtain error estimates that hold for arbitrary time displacements. 

The principal goal of this paper is to derive a modified version of Moyal quantum 
mechanics that treats the perturber degrees of freedom with a phase space formalism 
like (1.5)-(1.6) while maintaining a consistent 7ii-operator valued description of the 
radiator degrees of freedom. Within this hybrid operator-symbol formalism, the Moyal 
bracket induces a natural semiclassical expansion. The leading approximation is shown 
to recover the classical path approximation. Higher order corrections are well defined 
and sufficiently simple to allow numerical calculation. 



2. Mixed Weyl Symbol Calculus 

The mixed Weyl symbol is a parametric family of operators that arises when a phase 
space function is quantized in a subset of its canonical variables. The conventional 
symbol is a complex valued function in phase space, whereas the mixed Weyl symbol 
is operator valued. In this section, we modify the Wigner transform and Weyl 
quantization |]18[ procedures so that they define a mixed Weyl symbol formalism. Within 
this formalism we construct representations of the noncommutative star product, the 
Moyal bracket and the trace formulas that determine expectation values. The 
asymptotic expansions, which form the basis of a semiclassical analysis, are also 
obtained. 

Each of the Hi subspace components of Ti = Hi <8> r H- 2 has its own set of canonical 
operators and coordinates. We distinguish these variables by employing upper case 
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letters for Hi and lower case for H 2 . The d\ system classical phase space is T x * = 
T*(R§) ~ Rj x Rp 1 with coordinates Z = {Z x , Z 2dl ) = {Q x , . . . , Q dl ; P x , . . . , P dl ). 
In the Hilbert space Hi over R rfl , the quantized coordinate operators are Z = 
(Q,P) = {Qx, ■ ■ ■ , Qdi! Px, • • • > Pdi)- Acting on Rq 1 space wave functions $ G Hi, 
the Qj are the operators of multiplication by Qj and the conjugate momentum are 
Pj = —ihd/dQj. Likewise the d 2 system has phase space T 2 * = T*(R^ 2 ) with variables 
z = {qx, qd 2 ;Px, ■ ■ -,Pd 2 ) and canonical operators z = (q x , . . . , q d2 ;px, - - - ,Pd 2 )- 
The commutation relations 



[Z a , Zp 



ihJ (1) 



\Z a , z^ 







(2.1a) 



state the separate canonical character and mutual independence of Z and z. The 
matrices and are the standard symplectic matrices that arise on T x * and T 2 *. In 
block form 



J 



(i) 





-5° 



S di 




J(2) 







S d 2 





(2.16) 



where 5 di are the <ij-dimensional identity matrices. The transformations J are invertible 
with J' 1 = J T = — J. 



2.1. Partial Quantization 

In order to help formulate a statement of partial quantization, we recall the definition of 
the Weyl symbol appropriate for the full phase space T\* +2 = T*(Rq x R^ 2 ) ~ T x * x T 2 *. 



All the phase spaces Tj*,T 2 * and T* +2 are Euclidean. In this circumstance, the 
conventional Wigner transform [13, 19, || is well defined. It maps an operator A on H 
into a complex valued function that is supported on T* +2 . Specifically, if (X, x\A\X', x') 
denotes the coordinate space Dirac kernel of A, then the Weyl symbol is 

A W (Z; z)= [fdXdx e -i(r-x+P-*)/K(Q + iX,q + \x\A\Q -\X,q- \x) . (2.2a) 



Formula fl2.2ap is a restatement of ( |1.4[ ) adjusted to the larger phase space T* +2 . This 
(full) Wigner transform mapping A *—>■ A w is always denoted by a, independent of which 
Hilbert space H, Hi or H 2 that A acts on, i.e. A w = a A. 



The invertibility of the Fourier transform in (|2.2c| ) means that a has an inverse a 1 



The map o-~ l A v = A is Weyl quantization. The action of er -1 on exponential functions 
is a~ l [exp(iU ■ Z + iu- z)} = exp(iU ■ Z + iu-z), and may be viewed as the basis of Weyl 
quantization. Superposition of this result provides the quantization for arbitrary phase 
space functions. Suppose a is the Fourier transform of A w : T£ +2 —> C, then 



A W {Z; z)= J J dUdu a{U; u)e l(u - z+u - z) , (2.26) 
A = ffdUdu a{U; u)e i{L . (2.2c) 
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To verify that Weyl quantization ( 2.2<j) is the inverse of the Wigner transform ( p.2a| ), it 



suffices to take the \X, x) Dirac matrix elements of (|2.2d ). Fourier transform identities 
lead one to recover the Wigner transform (|2.2q ). 



It is evident from the structure of ( |2.2c| ) and the independence of the operators Z, z 
that one may partially quantize A w . In this context the two relevant choices are 

A wl (Z) = (a 2 l A^)(Z) = JJdU du a(U; u)e iU - z e iu - 7 , (2.3a) 
A w2 (z) = (a^A w )(z) = [fdUdu a(U;u)e iU Z e iu - z . (2.3b) 



In the notation above, the maps a 2 l and a^ 1 are the partial Weyl quantizations of A w 
with respect to the z and Z variables. The object A w i(Z) is an operator family on TC 2 
with a parametric dependence on Z. Likewise, A w2 (z) is an operator on TC\. 

The subscript labeling wl and w2 used in (2.3) is suggested by the Wigner transform 
point of view. For example, the wl subscript on A w i(Z) reminds one that the first 
argument of A, namely Z, has been dequantized to become a parametric dependence 
on Z. 

To pass from A wl or A w2 to A, one Weyl quantizes, respectively, either e lU ' z or e lu ' z 
in ( |2.1|) . After a Fourier transform this gives 

A = a^A wl = — jL_ JJdUdZ A wX (ZY u ^-^ , (2.4a) 
A = a^A w2 ee 77 ^— ffdudz A w2 (z)e lu ^- z) . (2.45) 



(2 



7T 



One may implement the transform A w t— > A via the a or b combinations of (2.3) and 
(2.4). This freedom to choose the order of the partial quantizations is a consequence of 
the mutual commutivity of Z and z. The product mappings that characterize these two 
equivalent quantization orderings are 

A = a~ 1 A w = cr ] ; 1 a 2 1 A w = a 2 l a^A^ . (2.5a) 

To complete the list of transformations we require the partial Wigner transforms 
<j\A = A w i and <j 2 A = A w2 . In view of ( |2.2q|) , (2.3) we have the definitions 

aiA = a 2 l aA, a 2 A = a^aA . (2.5b) 

It is straightforward to verify that aj 1 and <jj (j = 1, 2) are mutual inverses. The reason 
we have elected to define 0\ and a 2 by the product of transformations in ( |2.5fc| ) rather 



than by the direct route from A to A w % or A w2 via an integral like ( 2.2a ) is that the latter 
requires a non-standard form of the Dirac ket. Namely, instead of matrix elements of A 
with respect to the full ket \X, x) one would need matrix elements in terms of a partial 
ket | a;). 

The identification and interdependence of the various mixed symbols and partial 
quantizations are displayed in Fig. [TJ 
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It is instructive to contrast the basic features of the full Wigner-Weyl isomorphism 
a with those induced by the partial transformations <7j. Consider the tensor product 
operator A = f®g on T~C = Tii ® Ti.2- The o 2 transform of A has the factored form 
o~ 2 {f <8> g){z) = g(z)f, where g(z) = (ag)(z), cf. (1.3). In this example the z parametric 
dependence resides in the C- valued multiplier g(z). This means that the resultant 
families of operators commute, e.g. [g(z)f, g(z')f) = for all z,z'. However, for a 
general operator A, the different operators in the family A w2 (z) will not commute. 

A similar example concerns the a 2 image of / ® z a . Let n a be the linear coordinate 
functions on T 2 * , i.e. TT a {z) = z a . The 02 Wigner transform maps this operator into 

o- 2 {f®z a ){z) = (az a )(z)f = n a (z)f . (2.6) 

All the transformations a, o~i, 02 are linear bijective correspondences. A simplifying 
aspect of these transformations is the reality feature. Selfadjoint operators A = A< have 
real valued Weyl symbols A W (Z; z) = A W (Z; z)*, and Hermitian valued mixed symbols 
A w2 (z) = A w2 (zY. 

2.2. Symbol Products and the Moyal Bracket 

The * product of symbols [0, |2C], ^T| is constructed so that the noncommutative 
product structure of operators on Hilbert space is accurately mirrored in the space 
of symbols. For operators X, Y on 7i with symbols X w , Y w the conventional * product 
is X W *Y W = a{XY). By definition, this product is noncommutative. The mixed symbol 
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'star' is obtained via the same procedure. In terms of X w2 and Y w2 , let 

X w2 *Y w2 = a 2 (XY) . (2.7a) 

Of course, there is a different 'star' product for each of the symbols A w , A wl and A w2 . 
We denote all of these products by the same character *. The values of the surrounding 
symbols will determine the selection of the relevant 'star' operation. 



A basic representation of the U2-star product is given by the integral for inula ||22 



(\ \ 2(^2 
- J J J dz' dz" X w2 {z + z')Y w2 {z + z") exp [2i\(z' ■ J (2) z")] , (2.7b) 

where A = hT 1 . Whenever the integral exists, it provides an exact formula for X w2 *Y w2 . 
A distinguishing feature here is that the integrand is a product of two noncommuting 
Hi valued operators. This formula is verified by expressing XY in terms of two copies 



of (|2.4fc|) and using the Baker-Cambell-Hausdorff identity exp(iu' ■ z) exp(iu" ■ z) = 
exjp(—i~u'-J^u") exp(i(u' +u") -z) . Acting with a 2 on the XY product, while employing 
a 2 exp(i(u' + u") ■ z) = exp(i(u' + u") ■ z), leads to ( |2.7fi| ). Formula (|2.7?j ) also represents 
the a-star product if the quantities X w2 ,Y w2 are replaced with C- valued symbols. 

The mixed symbol Moyal bracket is the a 2 image of the commutator [X,y], 
specifically 

ih{X w2 , Y w2 } M = X w2 * Y w2 - Y w2 * A w2 = a 2 [X, Y\ . (2.7c) 

As with the * notation, the appropriate meaning of the bracket {•, -}m is selected by 
the values of its argument symbols. In all cases, the Moyal bracket is bilinear, skew and 
obeys the Jacobi identity. The Lie algebra defined by the Moyal bracket is isomorphic 
to the Lie algebra induced by the commutator of operators on 7i. 

2.3. Star Product Asymptotics 

In many circumstances the * product is close to commutative multiplication, and the 
difference of these two products is described by an asymptotic expansion involving 
derivatives of symbols. In this subsection we construct derivative expansions for mixed 
Weyl symbols and summarize their asymptotic structure. 

Asymptotic expansions of the * product are a consequence of the large A limit of the 
integral ( [2.7fc|) . In these expansions it is desirable to have an effective small h scaling, 
but our line shape application requires that Planck's constant be fixed at its physical 
value, h = 1.055 x 10 _34 J-s. In order to accommodate these opposing demands, we 
set A = (eh)^ 1 in ( |2.76|) , and effect the small Ti scaling by letting the dimensionless 



parameter e — > 0. 

First recall Groenewold's 121] expansion for C-valued symbols. Suppose the 
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operators /, g on TC 2 have symbols / = o~(f), g = o~(g). For small e, one has 

N-l I /j p fj\ n 

f * 9(z) = E - hr J ^ ■ ' ' J ^Lf^-M9^ n ^) + Rn{z) (2.8a) 

n=0 n - \ Z I 

Iff) 

= f(z)g(z) + —{f,g}(z) + 0(e 2 ). (2.8b) 

The first term on the right of (|2.86|) is ordinary multiplication, while the term linear in 
e is the Poisson bracket, {f,g} = V '/ '■ J^V g. The tensor indices [L\ ■ ■ ■ /J, n on / and g 
denote the partial derivative d n /dz^ 1 ■ ■ -dz^ n . The remainder Rn(z) is 0(e N ). 

We require expansions analogous to (2.8) in which / and g are replaced by mixed 
Weyl symbols. A useful derivative notation adapted to the Groenewold type expansion 
is the following. Denote by B the extended Poisson bracket operator for T 2 *. The B 
operator acts on a tuple of mixed Weyl symbols -< X w2 , Y w2 >~ to produce a new mixed 
Weyl symbol. The first and higher iterates of B are 

B ~< X w2 , Y w2 y (z) = J$X w2 ,^z)Y w2 , v (z) , (2.9a) 
B n -< X w2 , Y w2 y (z) = Jf) vi ■ ■ ■ J^%^ v ..MY^..^{z) (2.9b) 

= (-i) n 41 ■ ■ • • ( 2 - 9c ) 



The equivalent pair of representations ( 2.96|) and ( p.9cj) for B n follow from the skew 
nature of J( 2 \ namely Jj$ = —Jff. Clearly B -< X w2 ,Y w2 y (z) has the derivative 
structure of a Poisson bracket on T 2 , but with operator valued rather than scalar valued 
arguments. The quantity KB is dimensionless. 

The generalization of the Groenewold expansion to the mixed symbol * product is 

X w2 * Y w2 (z) = exp(iekB/2) -< X w2 , Y w2 y (z) (2.10a) 

= X w2 (z)Y w2 (z) + —B -< X w2 , Y w2 y (z) + • • • . (2.106) 

Formula (2.10) indicates how the * operator modifies the TCi product of two mixed 
symbols, X w2 (z) and Y w2 (z). Clearly, in the algebra of mixed symbols, X w2 (z), one 
has two noncommutative mechanisms — the * operation and the noncommutative 7ii 
operator product. 

As A —>■ 00 the exponential in ( 2.76 ) oscillates rapidly. This circumstance justifies an 



application of the stationary phase method |25j and leads to the asymptotic expansion 
( p,106|) . In order to establish the equivalent formal identity (|2 . 1 oj) , express the XY 



product as a multiple integral by using representation (|2.46|). In the integrand one 



encounters exp(— ^u- J^u') exp(i(u + u')-z). Rewrite this as exp(yV^ • J( 2 >V Z ') exp(iu- 
z + iu' ■ z')\ z=z i and take the left exponential operator outside the dudu' integration. 
This gives flnUfl) . 

An important application of (2.10) arises when the expansion is used to describe 
the o"2-Moyal bracket. Using the equivalence of fl2.96| ) and ( |2.9c| ) the exponential series 
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reorganizes as 

ih{X w2 ,Y w2 } M (z) = cos(yS) -< X w2 ,Y w2 y- (z) 

+ isiii(yB) -< X w2 ,f w2 y+ (z) (2.11a) 

= [X w2 (z),Y W 2(z)} + l -^B -< X w2 ,y w2 y + (z) + 0(e 2 ) . (2.116) 

Here -< X w2 , Y w2 >-±=-< X w2 , Y w2 y ± -< Y w2 , X w2 y denotes the symmetrized and anti- 
symmetrized tuple, respectively. Whenever the operator families X w2 (z) and Y w2 (z') 
commute, the cosine contribution vanishes and one recovers the familiar |]J sine function 
version of the Moyal bracket. 

The parameter e G [0, 1] measures the deformation of the * product from the 
conventional e = product. It provides a one dimensional variable that allows one 
to characterize the derivative expansions, (2.8) and (2.10), as e —>■ asymptotic 
approximations of the * operation. However, even if e = 1 and there is no exposed 
small parameter, the Groenewold formulas (2.8) and (2.10) may still be valid asymptotic 
expansions. (See Appendix A). 



2.4- Trace Identities 

Quantum expectation values are determined, within the density matrix formalism, by 
traces of operators on TC. For the a-transform the quantum trace is realized by a phase 
space integral. Assuming that Y and YX are trace class on H. with smooth symbols, 
the known f26| , Chap. IV] trace representations are 

T *nY = fi t dZ dz Y w (Z,z) , (2.12a) 

1+2 



h d ^+ d 2 JJt, 



Tr H YX=j^jJ T r/Z az YJZ. : )..Y W ( Z. : ) . (2/12/0 



1+2 



Our study of dynamics will systematically employ the mixed symbols Y w2 (z) and 
X w2 (z), thus we require the a 2 -counterparts of ( |2.4| ). One finds 



TthYX = -Lj^dz Tt Hi Y W 2(z) * X w2 (z) . (2.136) 

The proofs of (2.13) are simple. To show ( |2.13q| ), write the trace on Tii in terms of 
the Dirac integral kernel, viz. Tr^Y^z) = / dX (X\Y w2 (z)\X) . Now express Y w2 [z) 
via and use the fact that (X\ exp(z'C/'i ■ Q + iU 2 ■ P)\X) = exp(zC/a • X)5 dl {MJ 2 ). 

If one computes the right hand side of fl2.13ap the various Fourier integrals all become 
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delta functions, giving a result that is the dZ dz integral on the right of (|2.12a ). This 
establishes fl2.13a|) . 

Formula ( |2.13fc| ) is a consequence of ( p,13c| ), if X is replaced by YX in ( |2.13o| ). A 
further simplification of ( |2 . 1 3 allows one to replace the integrand by Tr^ i y w2 (^)A" w2 (z). 
In order to see this notice that an integration by parts shows 

J dz J^X^{z)Y^. v {z) = - J dz J$X w2 (z)Y w2 .^(z) = . (2.14) 

In the rightmost integral, the skew matrix Jj?J times the \iv symmetric operator Y w2 . >IJiV (z) 
sums to zero. A similar argument shows that the higher order B n terms vanish. 

The mathematical analysis of this section is heuristic and formal in style. Our 
intention has been to start from well known facts about the Wigner transforms and to 
keep the derivations at the simplest level consistent with obtaining all the structural 
formulas needed for the line shape problem. Nevertheless, we note that in the 
approximate theory given in Section 4, the space Tii is replaced by an iV-dimensional 
subspace 7i[ N ^ and so the mixed symbols become N x N matrix functions. For matrix 
valued symbols the rigorous methods and convergence estimates of pseudodifferential 
operator analysis (cf. Appendix A) will be applicable with straightforward modifications. 

A precursor of the mixed Weyl symbol formalism, employing matrix valued symbols, 
was used in |27j] to obtain equations of motion for quantum mean values in a Schrodinger 



evolution problem with a Hamiltonian containing spin structure. There one also finds 
an analog of the classical path approach. 



3. Interaction Picture Dynamics 

The final form of line shape theory we devise utilizes the interaction picture generated 
by the perturber Hamiltonian, H 2 . In this description an observable A acquires a 
time dependent form A(t) = T(H 2 ; r)A = U(H 2 ; t^AU(H 2 ; t). The operator T(H 2 ; r) 
denotes the H 2 Heisenberg evolution on Ti. When this interaction picture flow is stated 
in the o~ 2 -mixed symbol representation, it becomes 

a 2 (T(H 2 ;r)A)(z) = {r{H 2l r)a 2 A){z) = A w2 (r\z) . (3.1) 

Here r(H 2 ; r) = a 2 T(H 2 ; t)o~ 2 denotes the fundamental evolution operator that maps 
the initial mixed symbol A w2 (z) into its dynamical value A w2 (t\z). 

In this section, a semiclassical approximation of r(H 2 ;r) based on the quantum 
trajectories generated by H 2 = I\ <g> h 2 is constructed. The structure of the expansion 
closely parallels that found in Section I cf. (1.6) and in references 0, p^j . This prior work 
developed the semiclassical asymptotics for the (C-valued) Weyl symbol description of 
quantum mechanics. However, in the present circumstance A w2 (t\z) is an 7ii-operator 
valued symbol. The approach introduced in reference |3[ was to use quantum trajectories 
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as a basis from which to construct the full semiclassical expansion. The arguments below 
show that this idea remains applicable to the study of the mixed Weyl symbol evolution 
r(H 2 ;r). 

In Section 2, the dimensionless parameter e was introduced to describe the idea of 
scaling h to zero. We maintain this notation in this section and display this scaling as 
eh — > 0. The semiclassical representations found in this section presume that the a 2 
symbols for H 2 and A admit eh — > asymptotic expansions. This is indeed the case 
here, where both H 2 and A have h independent symbols. 



3.1. Quantum Trajectories 

Let the static canonical operators {z a }1 d2 be restricted to the space 7i 2 and let 
2q(t) = T(h 2 ;r)z a denote the associated h 2 Heisenberg evolution. The quantum 
trajectories are defined as z a (r\z) = r(h 2 ; r)ir a (z) = (z a (r)) w (z) and obey the equation 
of motion 

d 

—z a (r\z) = {z Q (r), h 2 } M {z) ■ (3.2a) 
at 

In this identity, {•, -}m is the Moyal bracket defined for C-valued Weyl symbols 
supported on T 2 *. At r = 0, the initial condition for z a is z a (0\z) = n a (z) = z a . 
Equation ( |3.2o| ) is an important example of ( |1.5fc| ). 

An approximate solution of z a (r\z) is obtained by expanding ( |3.2o| ) in small e. That 
this is possible is a result of the standard [0 small e asymptotic expansion of the Moyal 
bracket, cf. (2.8) 

00 (-l) n (eh /2) 2n 

K(r), h 2 } M (z) = K(r), h 2 }(z) + £ 1 I \ ' B 2n+1 -< z Q (r), h 2 y (z) . (3.2b) 

n=l \ An + l )- 

Using ( |3.2£|) in ( |3.2c| ) and solving for the coefficient functions order by order in e, one 
finds|| Sec - m ] 

*a{r\z) = t @£*?>{t\z) = g a (r\z) + ^z^\r\z) + 0(e 4 ) . (3.2c) 

n=0 U - Z 

The leading term above is the solution of Hamilton's equation 

g{r\z) = J^Vh 2 {g{r\z)) , (3.2(f) 

with initial condition g{0\z) = z. Because both the Moyal bracket and h 2 are even 
functions of h, one can prove that z a (r\z) (cf. || Lemma 5) is an fr-even function. As a 
result, all n odd terms vanish in expansion ( |3.2dJ . The higher order correction functions 
z^(t\z), n > 2, are obtained by solving an inhomogeneous Jacobi field equation, (cf. 
3.9^| ). The classical flow g(r\z) and coefficients z^(r\z) are defined and finite for all 
time displacements r. 
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3.2. Semiclassical Asymptotics of r(H 2 ;r) 

In this subsection, the eh — > asymptotic expansion for r(H 2 ; r) is derived. To start, 
we note that it is useful to reorganize the Fourier integral formulas for A and A w2 (z) in 
the following manner. Given A W (Z; z) with Fourier dual a(U; u), set 

A(u) = JdU a(U; u) e lU ^ . (3.3a) 

The quantity A{u) is a w-parameter function with TCi operator values. Representations 
( p.2cj) and ( |2.3t|) now read 



A = J duA(u)e luz , (3.36) 

A w2 {z) = J duA(u)e iu z . (3.3c) 
The action of T(H 2 ; r) on A has the integral form 

T{H 2 ; t)A = fdu A{u) e iu ^ r) . (3.3d) 



\z=0- 



Note that the integrand in ( ^.3d| ) may be restated as e tu ' z ^ = e ( 2; ( T )- s )- v ^ e m -( s + 2; )| 
Here s is an arbitrary vector in R^ 2 , and z(r) — s is shorthand for z(t) — sl 2 . Combining 
this with ( |3.3cj ) gives one an operator valued power series 



T(H 2 ; t)A = exp[(z(r) - s) ■ X7 z }A w2 (s + z 

OO 1 



z=0 

Yl f[( £ ( r ) ~ ' ' ' ( £ ( r ) ~~ s )Mn^w2; m ... M „(s) . (3.3e) 

71=0 U - 

Obviously, in this formula, A w2;(Ul ... Mn (s) is a static operator, while all the dynamics 
resides in the z(t) — s factors. Convert ( |3.3e| ) into an identity for r(H 2 ; r) by the 
application of the transform a 2 to obtain 

OO 1 

r(H 2 ; t)A w2 (z) = ^ -Mr) - s) m * • • • * (z(r) - s) Mn (z)l w2;w ... Mn ( S ) . (3.4) 

n=0 n - 



The series Q3.4|) is an s-dependent family of representations for the s-independent 
object r(H 2 ;r)A w2 (z). An optimal choice of s will make the nth series coefficient of 
order 0(e 1 ^), where l(ri) > n/2. For each given r, z choose s = z(r\z). This makes the 
n — 1 term zero, and gives us the expansion 

r(if 2 ; r)A w2 {z) = l w2 (z(r|^)) + Y, -W^ n (T\z)A w2] ^ n {z(r\z)) , (3.5a) 

n=2 n - 

= 5 n(z w (r) - s w ) * • • • * (z^(r) - s M J(z) | 8=b(tW . (3.56) 

Here 5 n is the permutation group operator (S^ = S n ) that averages over all the n! 
reorderings of the indices (/ii, • • • , « n ). The * products in ( |3.56| ) must be evaluated 
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before the constraint s = z(t\z) is imposed. In this evaluation the quantity s is a z- 
independent constant. The average S n ensures that the coefficient functions W are real 
and permutation invariant in the indices (/ii, • • • , 

An efficient way to compute the W functions is via the link operator L that expresses 
the extent to which the * product differs from commutative multiplication. This 
operator is Ljk = exp( 1 rr £?.,•&) ~ 1- It acts on an n-tuple of phase space functions 
-< A 1 {z^),A 2 {z^),---,A n (z^) y. The labels jk specify which pair of the n-tuple 
elements that B acts on prior to diagonal evaluation at 

z (l) = ... = z (n) = z _ 

The link 

operator is 0(e). A short calculation gives 
W/*w( r l*0 = L ^S2 -< z Ail (r),z M2 (r) y (z) 

= l -(eh) 2 B 2 12 S 2 -< z w (r),z w (r) y (z)+0(e 4 ) , (3.6a) 

lH2Ha( T \ Z ) — (-^12-^13-^23 ~ Li 2 L 23 )S 3 ~< Z Mi ( T ) ' Z M2 ( T ) 5 Z M3 ( T ) ^ ( Z ) 



1 



(efiyB 12 B 23 S 3 -< z„ 1 (r),z w (r) ) z w (r) 5- (z) + 0(e 4 ) . (3.66) 



A couple of remarks about representation (3.5)-(3.6) are in order. If A w2 (z) is 
a polynomial in z, then the series ( p.5a|) truncates at finite order and provides exact 



expressions for r(H 2 ;r)A w2 (z). In essence expansion (3.5) is a semiclassical expansion 
for an arbitrary mixed symbol A w2 . It employs z(r|z) (quantum flow) transport with 
higher order corrections arising from W tll ... tln {r\z). The link operator L was introduced 
in reference [[J in order to determine the symbol of an exponential operator exp A in 
terms of A w . 

3.3. The Standard Approximation 

Although expansion (3.5) is the basic semiclassical expansion for r(H 2 ;r), it is 
nevertheless not convenient for numerical calculation. If H 2 is z quadratic then 
z(t\z) = g(r\z) and W tM1 ... tMn {r\z) = 0. In this case the first term is exact, i.e. 
r{H 2 ;r)A w2 (z) = A w2 (g(r\z)). Generally, for n > 0, z^(t\z) and W w . 
not vanish and for increasing n these functions are difficult to numerically compute. 
In order to build a more computationally accessible approximation one combines 



expansion Q3.24 ) for z(t\z) with (3.5) and collects all terms of common e order. This 



results in an expansion of the form 

r(H 2 ; t)A w2 (z) = £ ^-L-^iH,; r)A w2 (z) . (3.7a) 

n=0 \ ln )- 

The coefficient operators ^ n \H 2 ] r) are composed of transport along g(r\z) plus a 
z-derivative structure inherited from ( |3.5dj . The first term is pure classical transport 

7 (0) (# 2 ; t)A w2 (z) = A w2 (g(r\z)) , (3.76) 
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while the 0(e 2 ) term reads 

^ 2 \H 2 - t)A w2 {z) = zW( T \z)A w2 .^g{T\z)) - ^w^{t\z)A w2 .^{ 9 {t\z)) (3.7c) 

+ j^w^ u x(r\z)A w2; ^x(g(r\z)) . 

The lower case w coefficients result from the suitably scaled limits of the functions 
W. Since W = 0(e 2 ) for Ws with two and three indices, one has 

w^(t\z) = \im(eh)- 2 W^{T\z) = B\ 2 -< g„(r),g u (r) y (z) , (3.8a) 



w 



t\z) = \im(eh) 2 W^ x (r\z) = B 12 B 23 -< g^r), g v (r), g x (r) y (z) 



- J {2) VgJr\z) ■ VVg u (r\z) ■ J^Vg x (r\z) . (3.86) 



These formulas make it evident that the w coefficients are functions of the classical flow 

g( r \ z )- 

A common building block in the expressions for ^ n \H 2 ; r) is the Jacobi field Vg(r\z) 
and its higher order derivatives Wg(r\z), etc. Jacobi fields describe the stability of a 
trajectory g(r\z) with respect to small changes of its initial data. Differentiating ( |3.2g| ) 
in the parameter z yields 



J{r)Vg{r\z) 



^-J (2) VVh 2 (g(r\z)) 



Vg{r\z) = . (3.9a) 



The solutions of ( |3.9g| ) with zero right-hand side are called Jacobi fields. The initial 
condition is Vg(0\z) = 5 (the 2d 2 identity matrix). Related functions Wg{r\z) solve 
a modified form of ( |3.9cj ) with a non-zero inhomogenous term |M Eq. (3.11)]. Finally, 
z( 2 \t\z) is the solution of the 2d 2 system of ODE's 

J{t)^z {2 \t\z) = - ^w^ 2 (T\z)J ( ^h 2mX {g(T\z)) 

+ Y^ w ^W2m( T \ z ) J !S h 2- m ^x(g{r\z)) , (3.9b) 

with initial conditions z^(0\z) = 0. 

We refer to expansion (3.7) as the standard semiclassical expansion of r(H 2 ; r). The 
basic structure of ^ n \H 2 ;r) is consistent with that derived in reference ||. The new 
feature here is that the target object A w2 is operator valued. In the case where the 
target symbol is C-valued, the existence of small e expansions of the standard form 
have been known for a considerable time. These alternate approaches [I5[ [16], [3], [L7[] 
are based on restructuring the Moyal equation of motion into a classical inhomogeneous 
Poisson bracket equation of motion. The inhomogenous component is formed from the 
non- leading terms of the derivative expansion of the bracket {•, -}m- Again one finds 
the leading 0(e°) term is classical transport. However, the higher order 0(e n ),n > 2, 
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terms have representations |3|, Eq.(4.16)] that are substantially more complicated than 
the ^ 2 \H 2 ;t) formula QE7jj ). 

Recently an extensive numerical study |28] of the effectiveness of the two term 
asymptotic expansion ( |1.66| ) was carried out. For systems with identical atom-atom 
pairs, such as helium, neon and argon, the time evolution of quantum expectation 
values was computed. The pair interaction in these systems was a phenomenologically 
determined Lennard- Jones potential. There was no difficulty in computing the functions 
(<7^(t|z), z^'(t\z), w^ u (t\z) and w^\(t\z)) which enabled the construction of ^°\h 2 ; r) 
and ^ 2 \h 2 \ r). The scattering problem was investigated for a variety of Gaussian 
initial states and observables. In most instances the x°\h 2 ; r) term dominated the 
contributions of ^ 2 \h 2l r) to the expectation value. 

One factor favoring the good convergence of expansion ( |1.6fc|) , in the problems of 
interest here, is that the classical system is completely integrable. This means that the 
classical flow is not chaotic. In particular, this implies that for almost all z G T 2 the 
Jacobi field Vg(r\z) can not have exponential growth in r. Nevertheless, there exists a 
set of unstable classical trajectories in this problem. All these motions are associated 
with an unstable equilibrium point. These points occur when the radial potential energy 
v e (r) = v(r) + L 2 /(2mr 2 ) (for a given angular momentum L) has v' e (r) = 0, with 
v"(r) < 0. Associated with this fixed point one obtains unstable orbits with positive 
Lyapunov exponent. The numerical studies we completed establish that the n = 2 term 
expansion ( |1.6£| ) is inaccurate for large time displacements in the small region of phase 
space surrounding the unstable fixed points . In this neighborhood of phase space, one 
must devise a different approximation for r(h 2 ;r). Several alternatives for overcoming 
this difficulty are offered in Section 5. 



4. Autocorrelation Representations 

The evolution T(H;t)p,j = U(H;typ,jU(H;t) controls the time behavior of the 
autocorrelation function, C(t). This section combines the mixed Weyl symbol formalism 
and the r(H 2 ; t) semiclassical expansion of Section 3 to construct numerically accessible 
approximate representations of C(t). 

First it is helpful to reformulate C(t), cf. (|1.3|), in terms of the H 2 picture evolution. 
Let U(t) = U (H ; t)U (H2; ty ; within the H 2 interaction picture framework, U(t) 
determines the full H dynamics. From the Schrodinger evolution equations for U(H; t) 
and U(H 2 ;t), it follows that U(t) is generated by the time dependent Hamiltonian, 

H{t) = T(H 2] -t){% + H 12 ) =% + T(H 2 ; -t)H X2 , (4.1a) 

ih^U(t) = U(t)H(t). {Alb) 

The operator U{t) is unitary and has initial condition £7(0) = I. The appearance of 
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H(t) to the right of U(t) is characteristic of a backward evolution equation. 

Conjugation of A by U(t) defines a Heisenberg evolution T(t) on TC in the standard 
way, T(t)A = U(t)^AU(t). Knowledge of the evolutions f(t) and T(H 2 ; —t) provides an 
alternate way of determining the correlation function. Employing the cyclic invariance 
of the trace and the commutation relation \p,j, H 2 ] = allows one to represent (|1.3| ) as 

C(t) = Tt h [(f (t)fij) (T(H 2 ; -t)p) ft] . (4.2) 

Whenever [p, H 2 ) = 0, the density matrix is time independent, i.e. T(H 2 ; —t)p = p. In 
most applications, the anisotropic interaction if 12 is small relative to H 1 and H 2 . This 
means that the radiator and perturber systems are weakly coupled and that the density 
matrix is accurately approximated by the tensor product of the form, p\{hi) <g> p 2 (h 2 ). 
For a system with radiators initially in the ground state, we may take the effective 
density matrix to be p = ®e _/3/l2 . We use this latter form of the density matrix 

throughout Section 4. 

4-1. Exact Mixed Symbol Dynamics 

The trace identity ( |2.13fi| ) determines C(t) in terms of the T 2 * phase space integral 
of (T(t)p,j) w2 (z) and (pp,j) w2 (z) = (e~ l3h2 ) w (z)\^ 1 ){(^ 1 \p, j . At this stage one needs 
an equation of motion for T(t) dynamics which is stated in terms of a mixed symbol 
evolution operator: (T(t)A) w2 (z) = r(t)A w2 (z), where r(t) = o 2 T{t)o 2 l . The symbol 
valued Heisenberg flow t{t)A w2 obeys 

-|r (t)A w2 = {T(t)A w2 , H(t) w2 } M (4.3a) 

with Moyal bracket (|2.7c|) . In the present application, A w2 is the dipole operator pj. 
The Hamiltonian here is the mixed symbol image of ( |4.1c| ), 

H(t) w2 (z) = H 1 + r(H 2 ; -t)(H 12 ) w2 (z) . (4.36) 

The matrix element version of (|4.3o|) with respect to the basis {1$^)}^=! reads 

a co 
ih-Xi(t\z) = J2{ X L(t\z) * ($ s \H(t) w2 (z)\$ r ) 

- ($ k \H(t) w2 (z)\<S> s ) * XMz)) . (4.3c) 

where X{ r {t\z) = (§ k \(r(t)pj)(z)\Q r ). This latter definition implies the ^-independent 
initial condition Xl r (0\z) = ($fc|/X,-| ( Iv). 

Given the solution of (^4. 3d), the correlation function is represented as 



C{t) = ^3 J dz Ti Hl (f{t)p,j){z) * {ptij) w2 (z) 
1 00 r 

= Ts E L dz XUt\z)Xl k (0\z)(e-^Uz) . (4.4) 



k=l 



r: 
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The second version of (fOj) uses the complete basis {l^)}^ to evaluate the Hi trace. 
Since the integrand is just the product of two symbols, cf. Q2.14|) , the * operation may 
be removed. Throughout this section the dimension d 2 = 3. 



4-2. Approximate Mixed Symbol Dynamics 



For typical molecular systems, equations ( |4.3c ) are too difficult to solve exactly. Further 
progress depends on developing approximate solutions that take advantage of the type 
of physics present in the line shape problem and the opportunities for simplification 
inherent in the mixed symbol formalism. The statement of exact dynamics provides 
a point of departure for the construction of various approximating methods. Possible 
options are the conversion of ( |4.3d.) into an integral equation, the use of perturbation 
theory, and the development of an eikonal representation. The approach we take here is 
to focus on the semiclassical structure. The subsequent reductions of (4.3) rest on three 
approximations: 

1) It is assumed that only a small number of h± eigenstates {\&k)}k=i significantly 
couple to each other. This means that 7i\ is replaced by the finite dimensional vector 



space, Hi"', spanned by the basis {\&k)}k=i- 



, opaiincu uy tuc u 

2) The Hamiltonian, H(t) W 2(z), may 
semiclassical expansion of Section 3, namely, 

oo ( eh y n 



be approximated by the standard 



H(t) w2 (z) = H 1 + J2 



-1 



(2n) 



n=o (2n)l 

Hi + (#l 2 )w2(<?Hk)) + 



(H 2 ;-t)(H 12 ) w2 (z) 
'eh) 2 



^ 2 \H 2 ;-t)(H 12 ) w2 (z) + 0(e 4 ) 



2! 



(4.5) 



3) The Moyal bracket on the right of fl4.3fl| ) can be replaced with leading terms 
of the small e expansion ( [2.116 ). 

Solvable reduced equations of motion are realized by applying approximations 1)- 
3). The finite coupled state assumption 1) is applicable when the range of thermal 
energies available in the collision process can excite a limited set of radiator eigenstates. 
The next stage is to employ the semiclassical expansion 2) valid to order 0(e 2 ). In 
this approximation, the H\ 2 part of H w2 (t) defines an 0{e l ) consistent time and state 
dependent molecular interaction by Mk r {t\z) = (&k\(Hi 2 ) w2 (g(— t\z))\$ r ). The modified 
version of ( [4.3 c]) thus becomes the N x N system 



i^xL(t\z) = (E r -E k )^ hr (t\z) 

N 



+ E {xis(t\z) * M sr (t\z) - M ks (t\z) * xiM< 



(4.6) 



The initial condition for 



is Xfcr(°l*) 



($*|&|$ r ). The (E r - E k ) X l(t\z) 
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contribution arises from the Hi part of H(t) w2 (z). The notation xKA z ) * s use d to 
distinguish the approximate solutions of ( f4.6|) from the exact X^(t\z). 

Observe that the non-zero term (E r — Ek)xi r {t\z) in ( |4.6|) means that the ih factor 
in front of the time derivative term can not cancel against a similar h multiplier on 
the right-hand side. In this sense ( |4.6| ) is very different in its h analytic structure than 
the Moyal equation of motion ( |1.5£j) , where such cancellation does occur. This means 
that methods for approximately solving ( |4.6| ) are not close parallels of the standard 
semiclassical expansions, ( |1.6fc| ) or ( |3.7o| ). 

The zeroth order reduction of (|4.6|) results if one replaces the * operation by ordinary 
multiplication to obtain the matrix ODE system 

l h^(t\z) = (E r -E k ) X ° k ^t\z) 

N 



+ E 



(4.7) 



[xli(t\z)M sr (t\z)-MMz)x°J(t\z) 

=1 L 

At this level of truncation one has recovered the classical path approximations to the 
line shape theory. The associated correlation function results from replacing X^(t\z) in 
( f4.4p with X°'''{'t\ z ) an d restricting the k sum to N terms. To further clarify this point 
consider again the rigid rotor example, cf. ( |1.2fc| ). Write the classical flow generated 
by h 2 in terms of its coordinate and momentum parts: g{t\z) = (q(t\z),p(t\z)). The 
molecular interaction M becomes 

oo a m=l 

M Ut\z) = E ^-T^lgHI^) £ Y^(q(-t\zm k \Y lm (Q)\^ r ) . (4.8) 

1=0 m=—l 

As is evident, the anisotropic interaction potential is evaluated along the classical 
path q(—t\z). The states {|$fc)} are eigenstates of the Hi system angular momenta, J 2 
and J z . This means that (&k\Yi m (Q)\& r ) has an explicit evaluation in terms of Clebsch- 
Gordan coefficients. Viewed as a phase space function M(t\z) is a time dependent 
valued symbol which is e independent. 

Now consider solutions of ( |4.6|) that include the leading noncommutative * effects. 
We organize this family of approximations as an asymptotic series in the small parameter 
e, which describes the deformation of the <T2-star product, cf. ( [2.10fc|) , about conventional 
Tii operator multiplication, 



Xl(t\z) = E <?J$(t\z) = x°d{t\z) + e X l d{t\z) + 0(6 2 ) . (4.9a) 



n=0 

The leading term Xk?(M z ) * s the classical path approximation. Placing ( |4.9 a|) in ( |4.6| ) 
and extracting the equation for x l, KA z ) gives 

<9 N r 

i^xltm = (E r - E k ) X %(t\z) + £ (xli(t\z)M sr (t\ z ) - M ks (t\z)xlr J (t\z)) 

01 s=l L 
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+ P^^/ k1m ^ Msr{tlz) ~ ^ a MkMz) ^ X -(t\z))} • (4.96) 

Here the initial condition is Xkr(^\ z ) = 0- The form of this equation for x 1,jf is a version 
of ( |4.7|) with an inhomogeneous term added. The J a l term on the right of ( 4.96 ) records, 
to leading order, the noncommutative nature of the * product. The z derivatives of the 
interaction M may be expressed via the Jacobi fields of the classical trajectories. By 
the chain rule one has 
d 



dz. 



-M sr (t\z) = ($ s \[V y (H 12 )„ 2 )](g(-t\z))\$ r )g 7 . a (-t\z) . (4.9c) 



The pattern one sees in (|4.96|) for the determination of x° S r an d Xsr continues 
to higher order. Given the values of {Xki}7=o the ODE system for Xkr results from 
combining fl4.6|) , ( |4 . 9 a|) with the * expansion (2.10). Normally, one would employ 
expansions 2) and 3) to a common order. Thus the 0(e 2 ) consistent calculation of xlr 
requires the addition of the j ( - 2 \H2; — t)(&k\(Hi2) W 2(z)\$ r ) contribution to Mk r (t\z). 

Evaluation of the integral ( |4.4j) is demanding because XkAA 2 ) needs to be 
numerically determined for each point z in the six dimensional phase space T 2 * . However 
the number of integration variables may be significantly reduced by using the spherical 
tensor structure present in this problem. The quantities (T(t)j2j) W 2(z) and (e~ l3h2 ) w (z)jlj 
are rank one tensors whose contraction is a scalar. One may reduce integration ( |4.4| ) 
to three parameters by representing z = (q,p) by the three Euler angles and three 
rotational invariants \q\,\p\ and q-p. The Euler angle integrals involve Wigner functions 
T> J M A 7) and can be done analytically. As an example of this type of phase space 
tensor reduction, see Section III]. 

The method of approximating X 3 kr (t\z) by the system of equations (4.9) depends 
upon representing the <T2-star product by its leading e terms. In the computation of 
observables one must have e = 1. As discussed in Appendix A this is an asymptotic 
derivative expansion. Roughly speaking, it will succeed if the higher order terms 
B n -< X J kr (t),M sr (t) y (z), n > 2 are small. 

It is useful to consolidate the correlation function results into a single statement. 
Assuming that xtr an d Xkr are solutions of ( f4.7|) and ( |4.9£|) , respectively, the 0(e x ) 
representation is 

I N r 

C(t) = T3 $>ilfcl**> / . dz (xli(t\z) + xHm)(e-^Uz) . (4.10) 

II k=l JT 2 

In our view the appropriate test of success of the approximating methods introduced 
here will be through numerical implementation and application to specific molecular 
systems. 

The correlation function formulas ( f4.4| ) and ( |4.10| ) assume that the perturber could 
be treated as a point particle. This restriction is easy to relax, by enlarging the 
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Hilbert space Ji\ to include both the radiator and perturber internal degrees of freedom. 



-0h2 



can take a variety 



Likewise, the presumed t = state of the system p 
of alternate forms which allow a superposition of molecular eigenstates {|$j)}] V and the 
replacement of e~ 2 by any function of h.2- The correlation function of interest here 
is based on the dipole operator fij. However, the method of this section continues to 
apply if p,j is replaced by any operator A having a fr-regularp4], |25|, [L6| mixed symbol. 



5. Conclusions 

The spectral profile of the intensity of an emission or absorption line is the Fourier 
transform of the dipole autocorrelation function Cjv(t). In this paper we have introduced 
a mixed Weyl symbol formalism to represent the dynamics needed to construct the 
single perturber C(t). Within this approach, the radiator-perturber relative separation 
variables are characterized by the phase space T 2 *, which serves as the support for the 
operator valued mixed symbols. Expectation values and in particular C (t) are realized, 



cf. (|4.4j) , by a phase space average over traces of the product of symbols. The result is a 
fully quantum theory of spectral line shapes. In the mixed symbol representation there 
is never any need to refer to perturber wave packets. 

Furthermore, the symbol equation of motion Q4.6| ) for the time evolving dipole 
moment, XkAA z )i admits a natural semi classical expansion which is based on the 
derivative expansion of the Moyal bracket. The mixed symbol formalism embeds within 
itself the classical path approximation, which appears as the leading order of the 
semiclassical truncation. The subsequent corrections arise from the noncommutative 
nature of the * product for the mixed symbol. 

The simplifications that result from the finite state coupling approximation 1) and 
the semiclassical expansions 2) and 3), give equations of motion that are suitable for 



numerical solution. To order 0{e 1 ) the equations ( 4/7) and (4.9) are no more elaborate 



than those previously used to numerically compute the atom-atom collision problem 
within the Moyal formalism, cf. [pgfl . The first term of fl4.10|) is entirely equivalent to 



the dipole autocorrelation function written in the classical path approximation and thus 
contains the same information as, for example, equation (4.2) in Griem's treatise |29f . 



It must be emphasized, however, that (|4.10|) gives the first two terms of a completely 
quantum mechanical expression for C(t). The established success of the classical path 
approximation in accounting for observed line shapes means that the 0(e°) version of 
the theory has much of the correct physics built into it. 

For heavy perturbers the classical path approximation works well. The diatomic 
molecular radiator-atomic perturber case has received much attention; for example, 



HCl-Ar collisions were studied by Nielsen and Gordon p0|. More recently, Looney and 



Herman pnj made a comparison of comprehensive calculations of the N 2 -broadened 
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rotational lines in the fundamental band of HC1 with state-of-the-art experimental data 
32] and found excellent agreement. Evidence that the 0{e 1 ) correction may be sufficient 



in many cases for an accurate quantum mechanical treatment comes from the work of 
Smith, Giraud and Cooper ][33| on CO-He. Their classical path calculations for this 
light perturber system agreed with close-coupling calculations to within 10%. 

Transforming the general autocorrelation function in ( |4.10| ) through ( |1.1| ) to give 
the profile I(oS) can be difficult. The line shape problem has two frequently employed 
limits, the impact approximation which applies to line centers and the quasi-static 
approximation which describes the far spectral wings |J. Various unified theories have 
been devised to connect these two regimes. Similar strategies will have to be applied to 
( |4.10| ) and (|1.1|) in order to compute actual profiles accurately. 

The mixed symbol method of computing C(t) has an important flexibility. The 
total isotropic intermolecular interaction energy is t>2 + V$ , i.e. the central potential 
can be arbitrarily divided between t> 2 and V-f^ . It is customary, in the classical 
path approximation, to set = 0. This means that the trajectories generated 

by h>2 = ft<2,o + ^2 are consistent with the full central potential. However, one can 
make other choices. If v 2 = 0, then h% = 7i 2) o (kinetic energy) and the classical 
trajectories are constant velocity straight lines. Here all the isotropic potential resides 
in V12 . For ft, 2 ,o evolution, the 0(e°) formula (4.5) in 2) is exact, namely H(t) w2 (z) = 
Hi + (Hi2) W 2{g{— t\z)). In this circumstance all the semiclassical behavior comes from 
approximation 3). This approach would improve upon the usual straight line trajectory 
technique often used to treat the broadening of atomic lines by ionic perturbers . 
Another appealing option is to place all the attractive part of the intermolecular 
interaction in Vj 2 ■ Now i> 2 w ih be purely repulsive and for this reason the h 2 flow 
will not have any unstable fixed points. In this way the unstable fixed point long time 
breakdown |28|, |34] , 17] of the standard semiclassical approximation is avoided. 

We have introduced the mixed Weyl symbol formalism in order to obtain a new 
and computationally viable full quantum version of line shape theory. However the 
mixed symbol representation is also applicable to any composite system having distinct 
quantum and semiclassical degrees of freedom. The scattering of spin dependent 
particles and the Coulomb excitation of an atom or a nucleus are additional examples 
of such systems. 
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Appendix A. Derivative Expansions 



The theory of pseudodifferential operators [go], |25| provides an analysis of the * product 
and its associated derivative expansions. For C- valued symbols /, g on T 2 * the * product 
is defined by the scalar analog of fl2.76p , namely 

/ * g = (nh)- 2d2 J J dz' dz" f(z + z')g(z + z") exp[2t(z' ■ J®z") /%] . (A.l) 

In the rigorous approach there are two important questions: a) The convergence of the 
integral ( |A.1|) and the status of the Groenewold series (2.8) as an asymptotic expansion 



of f*g; b) The properties of the operator corresponding to a given symbol, in particular 
the symbol which is the result of the * product and its approximations. 

The answer to these questions is effected by restricting the symbols to special classes 
of functions. A standard and important example is the class S m . Take (x,£) to be 
dimensionless versions of (q,p), and denote by S m a set of C°°(R d2 x TV 1 ' 2 ) functions 
with estimate 



V;!Vj/(.r.O < C ifc <0 mHil , (0 = V 1 + Kl 2 ( A - 2 ) 

valid for all multi-indices j = (ji, . . . ,jd 2 ), k = (ki, . . . , kd 2 ) and some constant Cjk, 
uniformly for x G H d ' 2 . 

Then the following statements 0, Theorem (2.49)] hold: 1) The map (/, g) — > / * g 
is continuous from S mi x S m2 -> S mi+m2 . 2) For / G S mi and g G S m2 } the remainder 
Rj\i in the Groenewold series belongs to the class S mi+m2 ~ N and so the expansion ( |2 . 8 oj ) 
becomes asymptotic. In other words, Rn = 0((^) mi+m2 ~ N ) and the error term vanishes 
as either N — > oo while (^) > 1, or as (^) — >• oo while N > mi + m 2 . In this sense 
the derivative series ( |2.8q|) is a valid asymptotic expansion even when the small scale 
parameter e is fixed at unity. 
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